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I .   INTRODUCTION 

In  the  opinion  of  the  author,  in  order  for  a  student 
to  realize  full  benefit  from  any  classroom  learning  process, 
that  process  should  contain,  whenever  possible,  some  type 
of  functional  visual  aid  or  laboratory  to  supplement  and 
support  the  in-class  lecture  periods.   The  purpose  of  the 
design  and  construction  of  the  stability  and  control  wind 
tunnel  model  described  in  this  thesis  was  to  provide  that 
additional  support  element  to  the  stability  and  control 
courses  offered  at  the  Naval  Postgraduate  School,  Monterey, 
California. 

Although  it  was  not  possible  to  incorporate  all  of  the 
degrees  of  freedom  of  an  actual  aircraft  into  the  model , 
the  model  was  designed  to  demonstrate  many  of  the  major 
dynamic  characteristics  of  a  full  scale  aircraft.   The 
response  of  the  model  differs  somewhat  from  that  of  an 
actual  aircraft,  due  to  the  fact  that  it  must  be  supported 
within  the  wind  tunnel  test  section  and  prohibited  from 
translating  in  any  of  the  three  coordinate  functions 
(Fig.  1-1).   This  difference,  however,  is  slight  and  does 
not  detract  from  the  primary  purpose  of  the  model. 
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Figure  1-1. 
Model  Reference  Stability  Axis  System 
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1 1 .   MODEL  DESCRIPTION 

A.   GENERAL  FEATURES 

The  model  (Fig.  1 1  - 1)  was  designed  to  illustrate  the 
characteristics  of  a  typical  single  engine  mid-wing  aircraft. 
With  minor  modifications,  however,  the  effects  of  various 
wing  geometries  can  be  studied  because  the  model  was  designed 
with  a  removable  wing  section.   Similarly,  model  tail  sizes 
and  shapes  are  variable.   By  the  use  of  various  wing  and 
tail  geometries  the  stability  characteristics  of  the  model 
may  be  changed. 

Model  attitude  was  controlled  by  the  deflection  of 
ailerons,  and  by  horizontal  and  vertical  tail  deflections. 
The  horizontal  and  vertical  tails  were  designed  as  Mall 
flying"  surfaces  so  as  to  eliminate  the  need  for  separate 
elevator  and  rudder  sections.   The  control  surfaces  are 
controlled  from  outside  the  tunnel  by  means  of  a  remote 
control  transmitter.   A  receiver  mounted  within  the  model 
receives  the  signals  from  the  transmitter  and,  by  means  of 
mechanical  linkage,  a  servo  deflects  the  appropriate  surface 
for  control  (Fig.  II-2). 

The  resulting  model  motion  was  detected  by  three 
potentiometers  which  were  situated  so  as  to  determine  the 
angular  rotation  of  the  model  about  the  X,  Y,  and  Z 
stability  axes.   The  potentiometer  output  was  connected 
to  a  chart  recorder  which  supplied  a  graphical  history 
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of  the  model  response.   The  motion  sensor  circuitry  is 
shown  in  Figure  II- 3. 

B.   PHYSICAL  MODEL  DESCRIPTION 

1 .  Fuselage 

The  fuselage  has  an  overall  length  of  29  inches 
and  was  of  circular  cross  section.   It  has  a  maximum  out- 
side diameter  of  4-1/2  -inches  and  a  fineness  rate  of  6.44. 
The  fuselage  is  constructed  from  a  block  of  mahogany  which 
was  hollowed  to  give  a  wall  thickness  of  1/4  inch. 

2.  Wing 

The  wing  is  constructed  of  mahogany  and  has  an 
NACA  2415  airfoil  section,  6  degrees  of  dihedral,  and  no 
geometric  twist.   The  wing  has  a  span  of  30  inches,  a 
leading  edge  sweep  of  23.5  degrees,  and  a  taper  ratio  of  .3 
The  left  and  right  wing  panels  are  attached  to  the  fuselage 
by  means  of  a  carry  through  section  (Fig.  II-4)  which  also 
houses  the  main  support  bearing  (Fig.  II-5).   The  rectangu- 
lar holes  in  the  carry  through  section  were  cut  to  reduce 
the  weight  of  the  unit. 

Ailerons,  which  are  8  inches  long  and  1  inch  wide, 
were  cut  from  the  trailing  edge  of  each  of  the  wing  panels. 
The  total  aileron  area  is  13.33  per  cent  of  the  wing  area. 

3 .  Horizontal  and  Vertical  Tails 

Both  horizontal  and  vertical  tail  surfaces  are 
constructed  of  mahogany  and  have  an  NACA  0012  airfoil 
section.   The  tails  are  designed  as  slab  surfaces  so  as  to 
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eliminate  the  necessity  for  separate  elevator  and  rudder 
sections.   The  tail  surfaces  are  held  stationary  on  rotating 
axles  by  means  of  set  screws.   The  axles  are  in  turn  sup- 
ported by  the  ball  bearing  tail  unit  assembly  shown  in 
Figure  II-6. 

The  horizontal  tail  has  an  overall  span  of  12.8 
inches  with  a  leading  edge  sweep  of  25.2  degrees  and  a 
taper  ratio  of  .24.   The  horizontal  tail  area  is  26  per 
cent  of  the  wing  area. 

The  vertical  tail  height  (from  the  fuselage  center- 
line)  is  6-1/2  inches,  the  leading  edge  sweep  is  50  degrees, 
and  the  taper  ratio  is  .31.   A  vertical  tail  area  to  wing 
area  ratio  of  .22  is  used. 

C.   MAIN  SUPPORT  BEARING 

The  model  is  supported  by  the  gimbaled  arrangement 
shown  in  Figure  II-5.   This  unit  is  installed  so  that  the 
pivot  point  of  the  bearing  is  located  exactly  14-1/2  inches 
aft  of  the  nose.   This  point  serves  as  the  origin  for  the 
model  reference  axes  as  well  as  the  model  center  of  gravity. 
Potentiometers  incorporated  within  this  unit  detect  model 
development  in  roll  and  pitch  up  to  angles  of  25  degrees. 
These  potentiometers  are  used  as  voltage  dividers,  and  their 
output  voltage  is  connected  to  a  chart  recorder  which  pro- 
vides a  graphical  measurement  of  the  model  response. 
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D.  TAIL  UNIT  ASSEMBLY 

The  assembly  shown  in  Figure  1 1 -  6  serves  to  support  the 
horizontal  and  vertical  tail  shafts  which  rotate  the  sur- 
faces to  the  desired  deflection  angle.   The  unit  is  con- 
structed of  aluminum  and  contains  four  ball  bearings  (2 
horizontal  and  2   vertical)  on  which  the  tail  axles  are 
supported.   The  entire  assembly  is  attached  to  the  fuselage 
at  a  point  such  that  the  horizontal  axle  is  exactly  12 
inches  aft  of  the  model  center  of  gravity.   This  resulted 
in  a  horizontal  tail  volume  ratio,  V„,  of  .732,  and  a 
vertical  tail  volume  ratio,  Vv,  of  .08. 

E.  WIND  TUNNEL  MOUNT 

The  model  was  supported  in  the  wind  tunnel  test  section 
by  means  of  the  mounting  assembly  shown  in  Figure  II- 7. 
This  support  was  constructed  of  aluminum  and  was  attached 
to  the  roof  of  the  wind  tunnel  as  shown  in  Figure  II- 8. 
It  has  an  overall  length  of  12  inches  of  which  approximately 
six  inches  protruded  down  into  the  tunnel  flow.   This  portion 
of  the  support  was  shrouded  by  a  symmetrical  aerodynamic 
fairing. 

The  model  was  connected  to  the  mounting  unit  through  a 
7/16  inch  diameter  stainless  steel  tube.   The  shaft  of  the 
main  support  bearing  was  inserted  into  the  bottom  of  the 
tube  and  tightened  in  place  with  two  set  screws.   Collars, 
in  conjunction  with  thrust  bearings  located  within  the  main 
unit,  absorbed  any  lift,  drag,  and  side  forces  of  the 
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model  while  still  allowing  it  to  rotate  freely  about  all 
the  coordinate  axes. 

The  shaft  of  a  potentiometer  located  on  top  of  the 
mounting  unit  was  inserted  into  the  top  of  the  stainless 
steel  tube  and  was  held  tight  with  a  set  screw.   This 
potentiometer  sensed  yaw  angle  of  the  model,  and  a  yaw 
stopper  attached  to  the  steel  tube  limited  yaw  to  30  degrees 

F.   SUMMARY  OF  MODEL  DIMENSIONS 

Wing  : 

Area  (including  portion  covered  j 

by  fuselage)  120  in 

Span  30  in 

Chord-root   '  6.15  in. 

-tip  1.85  in 

-MAC  4.38  in 

Aspect  Ratio  7.5 

Airfoil  Section  NACA  2415 

Leading  Edge  Sweep  23.5° 

Dihedral  Angle  6° 

Horizontal  Tail: 

Area  (including  portion  covered  2 

by  fuselage)  31.1  in 

Span  12.8  in 

Chord-root  3.94  in 

-tip  .93  in 

Aspect  Ratio  5.27 

Airfoil  Section  NACA  0012 


Leading  Edge  Sweep  25.2° 

Incidence  (variable)  10° 

Tail  Volume  Ratio  .732 

Vertical  Tail: 

Area  (including  portion  covered  2 

by  fuselage)  26  in 

Height  (from  fuselage  centerline)  6.5  in 

Chord-root  6.5  in 

-tip  1.5  in 

Airfoil  Section  NACA  0012 

Leading  Edge  Sweep  50° 

Tail  Volume  Ratio  .08 

Fuselage : 

Total  Length  29  in 

.  2 

Maximum  Cross  Sectional  Area  16  m 

Fineness  Ratio  6.44 
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III.   ANALYTICAL  MODEL 

The  model  was  designed,  and  its  motion  predicted  by 
means  of  linear  equations  of  motion.   From  these  equations 
it  was  possible  to  predict  model  response  to  various  control 
surface  inputs.   Also,  once  these  equations  had  been  estab- 
lished, it  was  possible  to  theoretically  vary  model  charac- 
teristics and  obtain  information  on  how  this  new  configuration 
would  respond. 

A.   LONGITUDINAL  EQUATIONS  OF  MOTION 

The  nondimensional  longitudinal  equations  of  motion 
were  obtained  from  Ref.  1.   These  equations  were  developed 
using  the  small  disturbance  theory  for  an  aircraft  having 
four  degrees  of  longitudinal  freedom.   The  model,  however, 
was  suspended  in  the  wind  tunnel  so  that  translations  along 
the  X  and  Z  axes  were  denied,  thereby  allowing  the  X  and  Z 
force  equations  to  be  neglected  entirely.   The  elevator 
hinge  moment  equation  was  also  neglected,  as  control  dis- 
placement rather  than  control  force  was  considered  as  an 
input. 

Due  to  the  assumptions  of  small  angles  and  constant 
airstream  velocity  in  the  tunnel,  perturbations  of  the 
velocity  in  the  X  direction,  U  terms,  were  zero.   Also, 
the  model  was  situated  within  the  tunnel  so  that  the  flight 
path  angle,  0,  and  the  angle  of  attack,  a,  were  equal. 
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Due  to  the  fact  that  the  horizontal  tail  was  a  slab 
surface  and  there  was  no  separate  elevator  section,  the 
quantity  n  which  was  originally  defined  as  elevator  per- 
turbation angle  was  replaced  by  tail  angle  of  attack  a    . 
It  was  further  assumed  that  a   terms  were  negligible. 

Applying  the  above  assumptions  to  the  initial  equations 
of  motion  yielded  the  following  longitudinal  equations  of 
motion  for  the  model. 


i  *a  -  (C    +  C   )  a  -  C    a-C     a.=  0 

B     K   m*  m  J  m       m    t 

a     q        a       a 


III-l 


q  =  a        1 1 1-2 

These  equations,  hoxvever,  do  not  take  into  account  the 
friction  inherent  in  the  main  support  bearing  on  which  the 
model  pivots.  Therefore,  they  must  be  further  modified  to 
account  for  this  frictional  retarding  moment. 

The  frictional  moment  present  at  the  pivot  is  a  function 
of  model  supporting  force.   More  accurately: 


M 


frict.  "  Vl    (W  "  L>  ITI-3 


where:   u-.  =  coefficient  of  friction  (in.) 

W   =  model  weight  (lbs.) 
L   =  model  lift  (lbs.) 

Putting  this  equation  in  nondimensional  format  results 
in  the  following  expression: 
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C   =     1    .  _1  c   a  1 1 1  -  4 

f   p  V  be    c    a 


3 
where:   p  =  air  density  (slugs/ft  ) 

V  =  velocity  (ft/sec) 


The  first  term  on  the  right  hand  side  of  equation  1 1 1 - 4  is 
the  constant  break-out  moment  of  the  bearing  for  any  given 
airspeed;  therefore,  the  perturbation  of  the  frictional 
moment  coefficient  during  motion  was  given  by: 


"yi 
AC    =  — -  CT   a  III -5 

f    c    a 


and  the  revised  model  longitudinal  equations  of  motion  are: 

yl 
iDa  -  (C    +  C   )  a  -  (C    -  — -  CT  )  a  -  C    a„_  =  0   III-6 
m«     m_        m_.       L./      m_.    t 

q  =  a   III -7 


a    c    a 


The  basic  rules  for  the  formulation  of  state  equations 
from  ordinary  differential  equations  [Ref.  2]  were  followed 
in  order  to  obtain  the  state  variable  representation  of 
the  model  longitudinal  equations  of  motion: 
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a 


(C    +  C   )/iD   (C CT  )/i 

v  m  m  • '   B    m^   -    Ly 
q     a         <X   c     a 

1  0 


B 


a 


(C    )/!„ 
^  m   ■"  B 


a 


III-8 


B.   LATERAL  EQUATIONS  OF  MOTION 

As  with  the  longitudinal  case,  the  lateral  nondimensional 
equations  of  motion  were  obtained  from  Ref .  1  and  then  modi- 
fied to  fit  the  explicit  conditions  of  the  model  suspended 
within  the  wind  tunnel.   To  begin  with,  the  Y  force  equation 
was  disregarded  due  to  the  fact  that  the  model  was  restricted 
from  translating  along  that  axis.   Also,  the  rudder  and 
aileron  hinge  moment  equations  were  neglected  since  control 
displacement  rather  than  control  force  constituted  the  input. 
The  term  z,    in  the  original  set  of  equations  denoted  rudder 
angle  perturbation;  however,  this  was  replaced  by  6   since 
the  vertical  fin  was  designed  and  constructed  as  a  slab 
surface.   It  was  further  assumed  that  all  £  and  6p  terms 
were  zero . 

The  model  was  suspended  within  the  tunnel  so  that  the 
rate  of  change  of  the  bank  angle,  $ ,  and  roll  rate,  p,  were 
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equal.   Also,  the  rate  of  change  of  the  sideslip  angle,  3, 
was  the  negative  of  the  yaw  rate,  r. 

Applying  the  above  conditions  to  the  initial  equations 
of  motion  resulted  in  the  following  model  lateral  equations 
of  motion. 


-C   3  +  (i.D  -  C   )p  -  (iED  +  C   )r  -  C  £,  -    C    6   =  0    III-9 

3            p             r  £      <5p 

-C   3  -  (iFD  +  C   )p  +  (i  D  -  C   )r  -  C   £  -  C    5   =  0    111-10 

n3      E     np       C     nr  n^     n6   F 


•    /\ 


(J)  =  p    III-ll 
3  =-r   111-12 


However,  as  before,  these  equations  do  not  reflect  the 
frictional  moments  present  and  must  be  modified  to  do  so. 
Model  side  force  is  the  major  parameter  contributing  to  the 
retarding  frictional  moments  during  model  lateral  motion. 
Therefore : 

Lr  .  „   =  ynY  111-13 

frict.     1 


and 


N,  •  .   =  y.Y  111-14 

frict.     2 


where:   Lf  •     =  frictional  component  of  rolling  moment 

(m-lbs  .  J 

N^  .     =  frictional  component  of  yawing  moment 
trie l •    /-•   ti„  "\ 
(m-lbs  . ) 
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u,  2  =  coefficients  of  friction  (in.) 
Y    =  side  force  (lbs.) 

The  frictional  coefficients  u,  and  y_  should  be  equal,  but 
the  coefficient  for  the  yaw  axis  has  a  different  value  due 
to  the  contribution  of  the  yaw  potentiometer. 

Putting  these  equations  in  nondimensional  format  results 
in  the  following  expressions  for  the  respective  frictional 
moment  perturbations. 

yl 
AC-   =  —■  C   6  111-15 

lf   b   y^ 

and 

y2 
AC    =  ~   C   3  111-16 

n£   b   y3 

where :   b  =  wing  span 

Incorporating  these  contributions  to  the  rolling  and 
yawing  moments  into  the  previous  set  of  equations  results 
in  the  following  corrected  set  of  model  lateral  equations  of 
motion. 

-  (C   +  A   C   )3  +  (iAD  -  C   )p  -  (i  D  +  C   )r 
13       y6  P  r 

-Cn  K    -    C   6,.  =  0       111-17 
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y9  ^ 

(C   +  rr-   C   )3  -  (irjD  +  C   )p  +  (i.D  -  C   )r 

nQ    b   yrJ  E     n/r     C     n' 

p       '  p  p             r 


C   £  -  C   6C  =  0 

*       6F 


111-18 


•       /s, 

(f)  =  p 
3  =  -r 


111-19 
111-20 


Putting  these  equations  in  standard  state  variable 
format  produced  the  following  matrix  system. 
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IV.   STABILITY  DERIVATIVE  RELATIONSHIPS 

In  order  for  the  equations  of  motion  to  be  of  any  use 
in  the  design  of  the  model,  the  individual  stability  deriva 
tives  were  programmed  and  expressed  in  terms  of  the  model 
physical  characteristics.   The  following  sections  present 
the  relationships  used  to  compute  the  required  stability 
derivatives . 

A.   LONGITUDINAL  STABILITY  DERIVATIVES 

The  values  of  the  stability  derivatives  involved  in  the 
longitudinal  system  of  equations  were  determined  using  the 
relationships  presented  in  Ref.  1.   These  are: 


C   =  -a(h   -  h)  IV-1 

C    -  -2a  f-^)Vu  IV- 2 

m       t  -    H 
q       c 


C    =  C   P-  IV-3 

m*     ni   da 
a     q 


C    =  -a+  Vu  IV-4 

m       t   H 

at 


whe 


re  V„  is  the  horizontal  tail  volume  ratio  given  by: 


v  =  z  z 
H   cS 
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From  the  above  equations  it  is  evident  that  the  longi- 
tudinal characteristics  of  the  model  are  very  dependent 
upon  the  tail  length  (1  )  and  the  tail  size  (S  ) ,  and  that 
by  varying  the  tail  area,  the  longitudinal  stability  charac 
teristics  of  the  model  may  vary  significantly.   For  the 
model  under  study  in  this  paper: 


C    =  -1.144  C   =  -14.979 

m  m 

a  q 


C   =  -6.231  C    =  -2.9448 

m*  m 

a  at 


B.   LATERAL  STABILITY  DERIVATIVES 

Unlike  the  longitudinal  case,  the  expressions  for  the 
lateral  stability  derivatives  are  very  complex  and  unwieldy 
to  use  and  program.   The  equations  of  Ref.  3  were  used  to 
estimate  these  derivatives  and  are  given  in  Appendix  A. 
For  this  particular  model ,  the  computer  program  solution 
had  the  following  values: 


Cn   =  -.1254  C   =  .07379 

X3  n6 


C,   =  -.4381  C   =  -.00251 


C,   =  .0437  C   =  -.15384 

r  r 
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C,   =  .69934 


C   =  -  .001197 


05554 


C    =  -.23747 
6F 
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V.   EXPERIMENTAL  LONGITUDINAL  RESPONSE 

A.  EQUIPMENT  CALIBRATION 

With  a  4-volt  potential  across  the  sensor  system,  the 
pitch  and  roll  potentiometers  were  calibrated  using  a 
clinometer,  and  the  yaw  potentiometer  was  calibrated  using 
a  protractor  with  the  yaw  stopper  as  a  pointer.   The  result 
ing  calibration  curves  are  presented  in  this  section  as 
Figure  V-l.   The  relationships  between  control  surface 
deflections  and  transmitter  stick  angle  were  also  measured 
and  are  presented  in  Figure  V-2. 

B.  TEST  PROCEDURE 

The  sensor  system  was  energized,  after  which  the  wind 
tunnel  was  activated  and  set  at  the  desired  airspeed.   The 
model  was  then  trimmed  out  at  a  slight  nose-up  attitude. 
The  horizontal  tail  surface  was  deflected  and  the  model 
response  to  a  step  function  input  was  observed  and  recorded 

Test  runs  at  airspeeds  of  50.7,  58.7,  65.4,  and  71.7 
ft/sec  and  tail  angle  changes  of  3,  4,  5,  and  6  degrees 
(positive  tail  angle  is  with  the  leading  edge  up)  were 
accomplished.   The  graphical  results  of  the  model  response 
to  a  tail  deflection  of  6  degrees  are  presented  as  typical 
in  Figures  V-3  through  V-6. 

From  these  results,  it  was  possible  to  determine  the 
values  of  the  damping  ratio,  £ ,  and  the  undamped  natural 
frequency,  oo  ,  for  the  model  at  the  different  airspeeds. 
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Figures  V- 7  and  V-8  clearly  illustrate  that  as  the  tunnel 
speed  is  increased,  the  value  of  oo   remains  fairly  constant, 
while  the  value  of  5  decreases.   This  results  in  a  gradual 
decline  of  the  damping  factor  2cw   as  the  airspeed  is 
increased.   This  decline  in  the  damping  term  is  a  conse- 
quence of  the  increased  wing  lift  at  the  higher  airspeeds 
which  reduces  the  pivot  force  resulting  in  less  bearing 
friction. 

Equation  1 1 1  -  6  was  used  to  obtain  the  characteristic 
equation  for  the  model  and  the  following  expression  for 

a)   and  r   were  obtained  as  functions  of  the  friction 
n 

coefficient ,  u-,  . 


(C   -  -JL  r   ) 

v  m     -    L 


w  =\i 2U_9 2L_  v-1 

n   W       iB 


-  (C   +  C   ) 
m*    m 

r  =  <L—. £L_  1  / _^ v-2 

21B      1/.C    ^C, 

m     -         L 
a   c    a 


Figures  V-9  and  V-10  illustrate  these  relationships  for 
this  model  assuming  constant  y-  values. 

The  experimental  results  were  used  to  determine  the 
unknown  bearing  friction  coefficient,  y, .   It  was  found 
that  a  value  of  1.75  was  required  to  match  the  predicted 
and  experimental  values  of  w  .   Figure  V-8  shows  the 
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resulting  comparison  values  for  the  damping  ratio.   It  is 
evident  that  as  the  tunnel  speed  increases  the  experimental 
and  theoretical  values  show  fair  agreement. 

From  the  various  test  runs  completed  it  was  found  that 
tunnel  speeds  around  65  ft/sec  exhibited  the  best  charac- 
teristics for  both  visual  observation  and  graphical  data 
acquisition  of  model  response.   At  velocities  greater  than 
75  ft/sec  considerable  model  buffeting  occurs  which  tends 
to  overshadow  any  data  obtained.   Conversely,  at  tunnel 
speeds  under  50  ft/sec  the  model  response  is  greatly  in- 
hibited by  bearing  friction,  and  the  resulting  heavily 
damped  motion  rarely  has  more  than  one  oscillation. 

C.   CONCLUSIONS 

The  model  studied  has  demonstrated  that  the  basic 
concept  of  a  remotely  controlled  stability  and  control  wind 
tunnel  test  model  is  a  valid  one.   It  has  been  shown  that 
for  the  longitudinal  case,  the  damping  and  frequency  of  the 
model  may  be  predicted  quite  accurately  using  linear 
equations  of  motion. 
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Figure  V-3.   Experimental  Angle  of  Attack  vs.  Time 
VQ  =  50.7  ft/sec,  Aat  =  6° 
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Figure  V-4.   Experimental  Angle  of  Attack  vs.  Time 
VQ  =  58.7  ft/sec,  Aa   -  6° 
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VI.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

Due  to  insufficient  time,  an  extensive  study  of  the 
lateral  characteristics  of  the  model  was  not  possible. 
However,  it  must  be  pointed  out  that  the  model  exhibited 
very  unstable  handling  qualities  in  roll.   In  fact,  this 
instability  was  detrimental  to  the  point  where  the  model 
had  to  be  restrained  from  rolling  in  order  to  obtain  the 
degree  of  accuracy  required  to  satisfactorily  determine  the 
model  longitudinal  characteristics.   Further  study  in  this 
area  of  the  model's  response  may  reveal  the  reason  for,  as 
well  as  a  solution  for,  this  unpredicted  instability.   The 
model  lateral  response  should  be  carefully  examined  and 
compared  with  the  theoretical  values  obtained  from  the 
computer  program. 

The  fact  that  there  presently  exist  difficulties  in 
the  roll  characteristics  of  the  model  should  not  immediately 
bring  a  conclusion  that  the  basic  concept  of  a  three- 
degree-of-freedom  stability  and  control  model  is  impractical. 
Rather,  with  continued  testing  and  refinement  it  is  believed 
that  a  tool  such  as  this  will  become  an  invaluable  asset 
to  students  actively  engaged  in  the  study  of  aircraft 
stability  and  control. 

Additionally,  it  is  suggested  that  the  friction 
coefficient  of  the  main  support  bearing  and  yaw  potentiometer 
be  experimentally  calibrated  so  that  the  model  response  may 
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be  correlated  more  accurately.   In  fact,  a  total  redesign 
of  the  main  support  bearing  for  less  friction  would  greatly 
simplify  and  increase  the  accuracy  of  the  predicted  response 
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APPENDIX  A 
LATERAL  STABILITY  DERIVATIVE  EQUATIONS 

The  following  equations  for  the  lateral  stability 
derivatives  were  obtained  from  Ref.  3,  and,  unless  other- 
wise stated,  any  specific  figures  mentioned  in  this  appendix 
are  from  that  reference. 

A.   VARIATION  OF  SIDE  FORCE  COEFFICIENT  WITH  SIDESLIP  ANGLE  • 
C 


y 


6 


C    =  C    +  C    +  C  A-l 

y8   ye    ye    y& 
p    pw    pb    pv 

The  wing  contribution  to  this  derivative  may  be  estimated 
from: 


C    =  -.0001  |r|  57.3  A-2 


where  T    is  the  wing  goemetric  dihedral  in  degrees. 
The  body  contribution  may  be  found  by  using: 

S 

C    =  -2  K.  Co2-)  A- 3 

y6       i  'S 


where:   K.  =  a  wing-body  interference  factor  obtained 
1    from  Fig.  7.1 

S   =  the  cross  sectional  area  of  the  fuselage  at 
the  point  along  the  body  where  the  flow 
ceases  to  be  potential. 

S   =   the  aircraft  reference  area. 


51 


The  vertical  tail  contribution  may  be  estimated  using: 

CyR   =  "k  CL    (1  +  3FJ  "V  T-  A"4 

Py  ay 

where:   k    =  an  empirical  factor  defined  in  Fig.  7.3 


C 


T    =  vertical  tail  lift  curve  slope 


aV 


s  /s       z 

(1  +  fe  n,.  =  .724  +  3.06  T— - — V,  ,,  +  .4  -^  +  .009  A 
v     d3   V  1  +  cosAc/4      d 


Sv   =  vertical  tail  area 

S    =  reference  area 

A    =  aspect  ratio 

Z     =  vertical  distance  from  the  wing  root  quarter 
chord  to  the  fuselage  centerline  (positive 
for  wing  above  centerline) 

d    =  fuselage  diameter  at  wing  intersection 


B.   VARIATION  OF  ROLLING  MOMENT  COEFFICIENT  WITH  SIDESLIP 

ANGLE  -  C 

X3 


C1   =  C     +  C     +  C  A-5 

6      6WB      3H      3V 


The  wing-body  contribution  may  be  determined  from: 
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=57.3  {CT    [(p-^O 


'WB 


LWR  -CL  'Ac/2  -MA  *f    ^CL  >A 


^   ^  +  (r— )a!  +  n^   KM 


AC. 


+  -r6!  +  (ACiA  } 

3  w 


A-6 


where:   CT    =  wing-body  lift  coefficient 
LWB 


Cl 

f— ^ 

lCL  JAc/2 


wing  sweep  contribution  obtained  from 
Fig.  7.11 


KM   =  the  compressibility  correction  to  sweep 
MA    from  Fig.  7.12 

Kf   =  the  fuselage  correction  factor  from  Fig.  7.13 


ci 


contribution  from  aspect  ratio  obtained 
from  Fig.  7.14 


T    =  dihedral  angle 


=  dihedral  effect  of  the  wing  from  Fig.  7.15 


K 


M, 


compressibility  correction  to  dihedral 
from  Fig.  7.16 


AC 


-=-£  =  -.0005  A  C^)2 


=  wing  span 


d    =  diameter  of  fuselage  at  wing  intersection 

p   w 
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Z  was  previously  defined  for  equation  A-4. 
w  l 

The  horizontal  tail  contribution  is  given  as: 


S    =  ci     !a  bii  a-? 

3H     3HB  S  b 


where:   S^  =  horizontal  tail  area 

n 

b„  =  horizontal  tail  span 

C,     is  obtained  by  solving  equation  A-6  using  the 
3WR  horizontal  tail  characteristics  vice  the  wing 

The  vertical  tail  contribution  is  given  by: 

r    =  r  (  v  cos  a    -   v  sin  a)         «  8 

where:   C    was  previously  determined  using  equation  A-4 


a  =  angle  of  attack  of  the  wing 

Z   =  absolute  distance  from  fuselage  centerline 
to  vertical  tail  aerodynamic  center  (a.c.) 

1   =  absolute  distance  from  model  e.g.  to  vertical 
tail  a.c. 


C.   VARIATION  OF  YAWING  MOMENT  COEFFICIENT  WITH  SIDESLIP 

ANGLE  -  C 

n3 

C   =  C    +  C  A-9 

nR    nft     nft 
3     3B     3V 

The  body  contribution  to  equation  A-9  is  given  by: 
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B   1 

Cnfi   "  "57-3  KNKRl  ST  IT  A"10 

BB  1 


where:   K„  =  an  empirical  factor  for  body  and  wing-body 
effects  from  Fig.  7.19. 

KR   =  Reynolds  Number  factor  for  the  fuselage 
Kl    from  Fig.  7.20 

SR   =  body  side  area 

1R   =  overall  fuselage  length 
The  vertical  tail  contribution  is : 

r  =  -f     (  v  cos  a  +   v  sm  a)       A-11 

where  all  the  elements  have  been  previously  defined. 


D.   VARIATION  OF  ROLLING  MOMENT  COEFFICIENT  WITH  ROLL 

RATE  -  C 

P 

C1   =  C1     +  C-    +  C1  A-12 

P     PWB     PH     PV 

The  wing-body  contribution  to  the  above  equation  is: 

ec. 


'1     l  k  }    3 
pWB     pW 


c1   :  c    M)  £  a- 13 


ecl 


where:   (  i  p)  =  roll  damping  parameter  from  Fig.  8.1 
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k  =  the  ratio  of  the  average  wing  section  lift 

curve  slope  C,    to  2tt 

aW 


J 2~ 

6  =  'l-NT 


The  horizontal  tail  contribution  is 


ci     ■  -5  cci  )H  sr  ^2  A"14 

ph  p 


where:   (C.  )„  is  obtained  by  solving  equation  A-13  for  the 
p   horizontal  tail. 

The  vertical  contribution  is  estimated  using: 

Z   , 
C1    =  2  G-^r  C  A-15 

pv         Yh 

where  all  the  variables  have  been  previously  defined  in 
this  appendix. 

E.   VARIATION  OF  YAWING  MOMENT  COEFFICIENT  WITH  ROLL 

RATE  -  C 
n 
P 


C    =  C    +  C  A-16 

n     n      n 

P     PW     Py 


The  wing  contribution  is  estimated  using  the  following 
equation. 
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c 

n 

C    =  -C1   tan  a  -  [-C1  tan  a  -  (p— E)      C,]  A-17 

nPw      Pw  P  L   CL=0 

M 


where : 

C-.    is  found  from  equation  A- 13 
PW 

a   =  wing  angle  of  attack 
C,  =  wing  lift  coefficient 

Cn  A  +4cosA  ,.   AB+i(AB  +  cosA  ,.)tan2A  ,„   Cn^ 

LL   C=0     AJ5  +  4COSA  /4   A  +1(A  +cosA   jtan-A       <-L   C  =0 
L  2         c/4       c/4        L 

M  M=0 

A  =  aspect  ratio 


2 


B   =Y1"m2cos"Ac/4 

M  =  mach  number 

2 
tanA  /A        tan  A  .. 

Cn  A+6(A+cosAc/4)  (|  —^  ♦    12  C/4) 

r       P~>  =    _    ±.  A  - 1  Q 

lCT     Jr    -n  6  A+4cosA    ~ 

L      CT=0  c/4 

M=0 


The  vertical  tail  contribution  may  be  found  from: 

9  Z  cos  a  -  1  sin  a  -Z 

C    =  -  f  (1  cos  a  +  Z   sin  a)  (— , — - -)  C    A- 20 

n^      b  v  v         v      J    K  b  vq 

Py  3V 

where  all  the  elements  of  A- 20  have  been  previously  defined. 
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F.   VARIATION  OF  ROLLING  MOMENT  COEFFICIENT  WITH  YAW 
RATE  -  C 


1 

r 


Cl   =  Cl   +  Cl  A"21 

r     rW     rV 


The  wing  contribution  is  given  by: 


C.  ACX 

ci    ■  cl  (cr^c=o +  c-r^  A"22 

rw        ML 


where:   Cy  =  wing  lift  coefficient 


l  •    AC1-B2)        +  AB+2cosAc/4  tan2Ac/4 

Ll         X  +  2B(AB+2cosA  ,J    AB+4cosA  ,„      8       1     A- 23 

( r.     =       *• c/4' c/4 r r-, 

Ll   LL=0  A+2cosA  ..  tan  A  ,,            ^L   CT =0 

M  i  +   c/4      c/4  L 

11  x    A+4cosA  ,.      8                  M=0 

c/4 


(p — )      is  obtained  from  Fig.  9.1 

L   V0 

M=0 


AClr  _  1    As"Ac/4 

r      12  A+4cosA  .. 

c/4 


T  =  wing  dihedral  angle  in  radians. 


The  vertical  tail  contribution  is  found  from 


58 


C,    =  -y  (1  cos  a  +  Z  sin  a)  (Z  cos  a  -  1  sin  a)  Cv      A-24 
rV   b  ^V 

where  all  the  variables  have  been  previously  defined. 


G.   VARIATION  OF  YAWING  MOMENT  COEFFICIENT  WITH  YAW 

RATE  -  C 

nr 


Cn   =  Cn   +  Cn  A"25 

rW     rV 


The  wing  contribution  is  found  from: 


Cn   =  <TT3  CL2  +  Ccr3  CD  A"26 

rW    CL  Do     ° 


where 


C 

( — j)      is  found  from  Fig.  9.4 
C 


C 

(p — )   is  found  from  Fig.  9.5. 
Do 


C^.     is  the  model  zero-lift  drag  coefficient. 
o 

The  vertical  tail  contribution  may  be  estimated  from: 

C    =  —  (1  cos  a   +  Z  sin  a)2  C  A-27 

n  r    b2    v         v         y6 
rV  PV 


where  all  the  variables  have  been  previously  determined 
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H.   VARIATION  OF  ROLLING  MOMENT  COEFFICIENT  WITH  AILERON 
DEFLECTION  -  C. 

The  procedure  for  determining  this  coefficient  is  a 
rather  complex  step-by-step  process  highly  dependent  upon 
the  exact  type  of  aileron  employed.   Therefore,  it  is  sug- 
gested that  Ref.  3  be  consulted  for  the  determination  of 
this  derivative. 

I.   VARIATION  OF  YAWING  MOMENT  COEFFICIENT  WITH  AILERON 
DEFLECTION  -  C 

C    =  K  CT  C.  A-28 

n^      L   1^ 

where:   K   =  an  empirical  factor  obtained  from  Fig.  11.3 
CL   =  the  model  steady  state  lift  coefficient 

C,   =  the  derivative  determined  in  section  H 

J.   VARIATION  OF  ROLLING  MOMENT  COEFFICIENT  WITH  RUDDER 
DEFLECTION  -  C, 

6F 

Reference  1  was  used  to  determine  the  value  of  this 

derivative.   The  following  equation  was  used. 


SV  Zv 

lr.        L     b    b 


where:   C.    =  vertical  tail  lift  curve  slope 


Sv   =  vertical  tail  area 
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K.   VARIATION  OF  YAWING  MOMENT  COEFFICIENT  WITH  RUDDER 
DEFLECTION  -  C 

6F 

Reference  1  was  used  to  determine  the  value  of  this 
derivative  also.   The  following  equation  was  used: 


S   1 

C    =  CT    ^-  ~  A-30 

nr     L    S   b 

6F     aV 


where  all  the  variables  are  known. 
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3}:  jjc  sjt>;: 

**  MODEL  CHARACTERISTICS  ** 


TOTAL    MODEL    PARAMETERS: 


TRIM    LIFT    COEFFICIENT     (CD     = 


0.0 


MODEL    ANGLE    OF    ATTACK    (DEG.)     =    -5 . 0982669D-01 


TAIL    ANGLE    OF    ATTACK     (DEG)     = 


2.4765951D    00 


TOTAL    LENGTH    (IN)     = 


2.9000000D    01 


TOTAL    WEIGHT    (LBS)     = 


4.6997845D    00 


STICK    FIXED    STATIC    MARGIN    = 


2.2599734D-01 


LIFT    CURVE    SLOPE     (PER    RAD.)     =         5.062'+000D    00 


MOMENTS    OF     INERTIA     (IN2-LBS) 


I(X)     =       6.2023497D     01 


I (Y)     =  2.6436234D    02 


I(Z)     =       3.1432216D    02 


KXZ)    =        1.0348581D    01 


MODEL    C.G.    WITH    RESPECT    TO    REFERENCE    AXIS    ORIGIN    (IN) 


XBAR    =       4. 7068655D-02 


YBAR    =      -3.1916357D-03 


Z8AR    =    -K6414961D-02 
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WING  AND  TAIL  SURFACE  PARAMETERS: 

-WING  CHARACTERISTICS- 
ROOT  CHORD  (IN)  =  6.1500000D  00 
TIP  CHORD  (IN)  =  1.8500000D  00 
WING  SPAN  (IN)  =  3.0000000D  01 
WING  AREA  (IN2)  =  1.2000000D  02 
L.E.  SWEEP  (DEG)  =  2.3500000D  01 
DIHEDRAL  ANGLE  (DEG)  =  6.0000000D  00 
LIFT  CURVE  SLOPE  (PER  RAD)  =    4.5309000D  00 

-HORIZONTAL  TAIL  CHARACTERISTICS- 
ROOT  CHORD  (IN)  =  3.9375000D  00 
TIP  CHORD  (IN)  =  9.3000000D-01 
TAIL  SPAN  (IN)  =  1.3125000D  01 
TAIL  AREA  (IN2)  =  3.1100000D  01 
L.E.  SWEEP  (DEG)  =  2.5200000D  01 
DIHEDRAL  ANGLE  (DEG)  =  0.0 
LIFT  CURVE  SLOPE  (PER  RAD)  =    4.0211000D  00 
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-VERTICAL  TAIL  CHARACTERISTICS- 
ROOT  CHORD  (IN)  =  6.5000000D  00 
TIP  CHORD  (IN)  =  1.50000000  00 
TAIL  HEIGHT  (IN)  =  6.375000CD  00 
TAIL  AREA  ( I N2 )  =  2.6000000D  01 
L.E.  SWEEP  (OEG)  =  5.00000000  01 
LIFT  CURVE  SLOPE  (PER  RAD)  =    2.989100CD  00 
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LONGITUDINAL  STABILITY  DERIVATIVES: 

CM-  =  VARIATION  OF  PITCHING  MOMENT  COEFFICIENT 

WITH  - 
A  =  ANGLE  OF  ATTACK  (ALFA)     Q  =  PITCH  RATE 
A(DOT)  ---  VARIATION  OF  ALFA  WITH  RESPECT  TO  TIME 


CMQ  =       -1.4979070D  01     CMA  =        -1.1440915D  00 
CMA(DOT)  =  -6,2308383D  00     CMA(TAIL)  =  -2.9448839D  00 

ELEMENTS  OF  STATE  VARIABLE  MATRIX  A  ARE: 

-4. 49702340-03      -2 . 42575 59D-04 
l.OOOOOOOD  00       0.0 

ELEMENTS  OF  THE  B  VECTOR  ARE: 

-6.2438798D-04 
0.0 
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LATERAL    STABILITY    DERIVATIVES: 

CL-    =    VARIATION    OF    ROLLING    MOMENT    COEFFICIENT 

WITH    - 
CN-    =    VARIATION    OF    YAWING    MOMENT    COEFFICIENT 

WITH    - 
B    =    SIDESLIP    ANGLE     (BETA) 

P    =    ROLL    RATE  R    =    YAW    RATE 

DA    =    AILERON    DEFLECTION  DF    =    RUDDER    DEFLECTION 


CLB    =      -1.2537315D-01 


CNB    =  7.3790651D-02 


CLP    =       -4.3806752D-01 


CNP    =       -2.5131716D-03 


CLR    =  4.3663509D-02 


CNR    =       -1.5334192D-01 


CLDA    =       6.9936720D-01 


CNDA    =    -1.1977032D-03 


CLDF    =       5.5541463D-02 


CNDF    =    -2.3746739D-01 


ELEMENTS    OF    STATE    VARIABLE    MATRIX    A    ARE: 


1.2747884D-01 
4.3405394D-03 
0.0 
l.OOOOOOOD    00 


1. 1230146D-02 
-8.4139888D-03 
-l.OOOOOOOD    00 

0.0 


-3.5770197D-02  0.0 

3.0354552D-03  0.0 

0.0  0.0 

0.0  0.0 


ELEMENTS  OF  THE  B  MATRIX  ARE: 


2.0346785D-01 

6.6304872D-03 

0.0 

0.0 


1.3884961D-02 
-1.3101244D-02 
0.0 
0.0 
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